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We introduce an analytical approximation scheme to diagonalize parabolically confined two di- 
mensional electron systems with both the Rashba and Dresselhaus spin-orbit interactions. The 
starting point of our perturbative expansion is a zeroth-order Hamiltonian for an electron confined 
in a quantum wire with an effective spin-orbit induced magnetic field along the wire, obtained by 
properly rotating the usual spin-orbit Hamiltonian. We find that the spin-orbit-related transverse 
coupling terms can be recast into two parts W and V, which couple crossing and non-crossing ad- 
jacent transverse modes, respectively. Interestingly, the zeroth-order Hamiltonian together with W 
can be solved exactly, as it maps onto the Jaynes-Cummings model of quantum optics. We treat 
the V coupling by performing a Schrieffer- Wolff transformation. This allows us to obtain an effec- 
tive Hamiltonian to third order in the coupling strength of V, which can be straightforwardly 
diagonalized via an additional unitary transformation. We also apply our approach to other types 
of effective parabolic confinement, e.g., 2D electrons in a perpendicular magnetic field. To demon- 
strate the usefulness of our approximate eigensolutions, we obtain analytical expressions for the n th 
Landau-level (j„-factors in the presence of both Rashba and Dresselhaus couplings. For small values 
of the bulk g-factors, we find that spin-orbit effects cancel out entirely for particular values of the 
spin-orbit couplings. By solving simple transcendental equations we also obtain the band minima 
of a Rashba-coupled quantum wire as a function of an external magnetic field. These can be used 
to describe Shubnikov-de Haas oscillations. This procedure makes it easier to extract the strength 
of the spin-orbit interaction in these systems via proper fitting of the data. 

PACS numbers: 72.25.Dc,71.70.Ej,73.63.-b 



I. INTRODUCTION 

For a wide range of systems studied in quantum trans- 
port the starting point of the sample fabrication is a 
two dimensional electron gas (2DEG). Various nanostruc- 
tures are defined in the 2DEG by using either metallic 
gates to expel electrons or by etching into the electron 
gasi. Quantum wires can be formed using both methods 
and in some cases the wire confinement can be assumed 
parabolic. 

Since the seminal experiments showing conductance 
quantization—!. 3 in narrow gate-tunable constrictions act- 
ing as quantum wires (or quantum point contacts) de- 
fined in 2DEGs, these structures have been extensively 
used to investigate a rich variety of physical phenomena, 
e.g., the 0.7 anomaly and analogs^ 5 -. Spin orbit interac- 
tion has also been investigated experimentally-^ in par- 
allel quantum wires, where universal conductance fluctu- 
ations are suppressed. More recently, ballistic spin reso- 
nance due to an intrinsically oscillating spin-orbit field 
has been experimentally realized in a quantum wire 8 . 
The observation of the 'one-dimensional spin-orbit gap' 
in quantum wires has also been reported recently^. Both 
of these experiments highlight the use of spin-orbit ef- 
fects in quantum wires as a means to control the spin of 
carriers, an important ingredient for potential spintronic 
applications. 



Bulk semiconductors lacking an inversion center in the 
crystal lattice have a built-in spin-orbit interaction, the 
so called Dresselhaus termiS. Heterostructures in which 
the 2DEG is formed by an asymmetric confining po- 
tential also exhibit the Rashba spin-orbit interaction^ 1 -. 
As shown experimentally, the Rashba coupling strength 
can be tuned via proper gating of the structur e 12 ' 13 , 
which makes the Rashba interaction very appealing for 
potential technological applications involving spin con- 
trol. More recently, yet a new type of spin-orbit interac- 
tion has been found in symmetric two-dimensional quan- 
tum structures with two subbands: the inter-subband- 
induced spin-orbit couplin g 14 ' 15 . Interestingly, this new 
spin-orbit interaction gives rise to non-zero intrinsic spin 
Hall effect^. 

Here we consider a parabolically confined asymmet- 
ric (narrow) quantum well with only the lowest occupied 
subband and in the presence of both the Dresselhaus and 
the Rashba terms 

H = T^ipl + P 2 y) + J-(Py< J x ~ P* a y) + 

-(p x <T x - Py<7 y ) + ^m*ujly 2 , (1) 

where m* is the electron effective mass, ujq characterizes 
the strength of the parabolic confinement, p x t y \ is the 
momentum operator in the x(y) direction, and a and /3 
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are the Rashba and the Dresselhaus coupling strengths, 
respectively. Similar systems have been studied before 
in the presence of Landau level quantization, using ei- 
ther variational (Hartree-Fock) methods^!, second order 
perturbatio n 18 ' 19 or by obtaining the energy spectrum 
in terms of perturbation series^. However, we are un- 
aware of any analytical solution for both Rashba and 
Dresselhaus that goes beyond second order perturbation 
for multiple levels, i.e., not just treating the two low- 
est orbital levels as in Ref^.) We note that in the ab- 
sence of magnetic fields the spectrum of 2D electrons with 
both the (linear-in-momentum) Rashba and Dresselhaus 
interactions is known; for non-zero magnetic fields only 
particular cases are solvable but the general solution is 
not known. The above model Hamiltonian has also been 
investigated in connection with the Zitterbewegung of in- 
jected spin-polarized wave packets in quantum wells and 



wires 
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We have developed a perturbative scheme to diago- 
nalize our Hamiltonian Eq. (TTJ) based on separating it 
into terms that can be treated exactly ( "effective Jaynes- 
Cummings model" ) plus a perturbative term . In Section 
HIl we derive the perturbative scheme in the case of a 
parabolic quantum wire with Rashba spin-orbit coupling 
and in Scction lHTl we extend it to include both the Rashba 
and Dresselhaus interactions. These results are then ap- 
plied in Section IIVI to magnetotransport in (i) quantum 
wires with Rashba coupling and (ii) 2DEGs in a perpen- 
dicular magnetic field in the presence of both Rashba 
and Dresselhaus couplings. Both of these systems can be 
mapped onto the Hamiltonian in Eq. (1) and the general 
analytical results derived straightforwardly used. For (i) 
the approximate eigenvalues allow us to obtain beating 
patterns that directly relate to the Shubnikov-de Haas 
oscillations, by simply finding the roots of transcenden- 
tal equations. This can be used to more easily extract 
the Rasbha coupling strength from experimental data 6 . 
For (ii) we derive an expression for the effective (/-factor 
of Landau levels. Interestingly, we find that when the 
g-factor is small enough, the spin-orbit effects cancel out 
exactly for certain values of the Rashba and Dresselhaus 
couplings. 



II. PURE RASHBA SPIN-ORBIT COUPLING 

When only the Rashba term is present the Hamiltonian 
of the system is 

It is convenient to introduce the standard ladder operator 
which yields the new Hamiltonian 



H n = 



Pi , t t hk nPx 
+ hLUoa' a ov 



2m 



ih 2 kji 



where we have introduced the Rashba wave vector kR = 
^p 2 , the oscillator length £ = \J h/m*uiQ, and the energy 



is measured relative to fuvo/2. The system is translation- 
ally invariant along the wire (i.e., \p x ,Hn] = 0) so we 
seek eigensolutions of Hr which are plane waves in this 
direction and hence eigenvectors of p x with eigenvalues 



Hr = -k 2 + a 1 a 



k R ka y + "y=( at ~ a ) a x- ( 4 ) 



The energy in Eq. Q is measured in units of hujQ, and 
momenta (both hk and hk^) in l^ 1 . At a first glance, the 
above Hamiltonian looks like a shifted harmonic oscilla- 
tor; however, since a x and a y do not commute, an exact 
solution is not known. 

Our goal is to construct a perturbation expansion in 
the small parameter kR but note that the product k^k 
need not be small. It is convenient to rotate the spin 
operators to obtain a new Hamiltonian where the effec- 
tive magnetic field due to the momentum along the wire 
couples to a z 



= — k 2 + a^a — /cR,fcer z + — p(<7+(a^ — a) — h.c. ), 
2 v 8 

(5) 

where er + = a x + ia y . Important features of the spectrum 
and the strategy of the following perturbative expansion 
are best explained in terms of the zeroth-order Hamilto- 
nian 



Hq = — k 2 + a)a — kRka z . 



(6) 



The kets \k,n,s) are eigenstates of H . These are also 
simultaneous eigenstates of a^a and a z with eigenvalues 
n = 0,1,... and s — ±1, respectively. The eigenval- 
ues £° s (fc) — \k 2 = n — skRk are plotted in Fig. Q] for 
Jcr = 0.25. Focusing on positive wavevectors k > one 
sees that the state \k,n, —1) crosses the state \k, n+1, +1) 
at k* — This means that non-degenerate pertur- 

bation theory cannot be applied to the terms contain- 
ing a)<j+ (or their Hermitian conjugate) as they couple 
these states. On the other hand, states \k,n, +1) never 
cross \k,n+ 1,-1) and thus the coupling terms aV_, 
and their hermitian conjugate, can be treated perturba- 
tively. The negative wavevector part of the spectrum is 
analyzed in the same manner using the Kramers relation 



-k) 



,(k), i.e., for k < the spin indices in the 



above discussion are reversed. 

The above considerations motivate us to separate the 
transverse mode coupling into two parts W and V 



W 



V 



— (a'cr+ - ao-), 
ikR, 



V8 



— acr + ), 



(7) 
(8) 



where W contains the term that couple states that cross 
and V is the part of the coupling that can be treated 
pertubatively (i.e., V does not couple states that cross). 
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where S is chosen such that 




FIG. 1. Energy spectrum of the zeroth-order Hamiltonian 
Ho-\k 2 [Eq. ©] for fc R = 0.25 (first ten wire modes shown). 
The circles and squares indicate the crossings between adja- 
cent and next-to-adjacent, respectively, opposite-spin oscilla- 
tor states. Here k* = defines the crossing point of the 
energy dispersions of the states \k, n, —1) and \k, n + 1, +1). 



It is worth pointing out that V and W are related via 
time reversal symmetry, i.e., V = T~ 1 WT 1 where T is 
the time reversal operator—, i.e. their roles are reversed 
for the negative k states. The following approximation 
scheme docs not rest upon time reversal symmetry. Time 
reversal is only used to obtain the negative wavevector 
(k < 0) states from the positive ones. 

Thus we split the Hamiltonian in Eq. ([5]) into two 
parts: 



Hti — H] 



R,JC 



V, 



(9) 



where V is defined in Eq. ((5J) and H R jc is a Jaynes- 
Cummings-like Hamiltonian—^ 



H n , jc = H + W 

1 „ . 

= — k + a' a — k R ka z 



ik R , t , 
-j=(a*(j + - acr_), 



(10) 



which is exactly diagonalizable. In fact, we show in 
Appendix [X] that a generalized coupling of the form 
[ r y(n)cL ! (T+ + h.c], with 7 being some generic complex 
function of h = a^a, can be diagonalized. This will guide 
us in determining which terms to keep and which to dis- 
card in the following perturbative expansion in V. 



A. Effective Rashba wire Hamiltonian 

We start with defining a new Hamiltonian using the 
transformation 



H R = e s H R e- s , 



(11) 



[S,H + W] + V = + O(k R ). 



(12) 



This procedure is the same as the usual Schrieffcr- Wolff 
transformation—, the only difference is that the zeroth- 
order Hamiltonian in Eq. (|10[) contains a term propor- 
tional to the perturbative parameter k R . 28 . We find that 



S 



-ik R a)a 



V8(l + 2fc R fc) 4(1 + 2k R k) 



k" R Akaa + a i a 



kjB k 



2k R k 



-a' a a 



h.c. (13) 



In the last line we have introduced the factors Ak and 
Bk which depend only on k. Here we do not write out 
the explicit forms of Ak and Bk since they will only oc- 
cur in terms in the effective Hamiltoninan oc k R (the 
extra power of k R coming from [S, V]) and can thus be 
discarded when certain restrictions on k are taken into 
account. The denominator in the £?fc-term vanishes when 
k = 3k*, since at that k- value it couples states that are 
degenerate. In order to treat this term perturbatively we 
make the restriction |3fc* — k\ > k R . 

We can then proceed with the usual calculations, i.e. 
knowing the form of S the transformed Hamiltonian is 



H R = H + W 
+0(4). 



1[S,V] + ±[S, [S,V]} 



(14) 



Inserting Eq. (p~3|) into Eq. (fl4|) we obtain the following 
exact result (up to the order indicated) for the Hamilto- 
nian 



1 - T z (2a' f a + 1) 
Mi ~"" " " l! 4(l + 2fc R fc) 



HL= H + W + kl 



+k- 



K 8V2{l + 2k R k) 
1 



{a/ aa) a + — h.c.) 



+k R C k 



t 1+ 3 

aa+a' a -\ — a 1 cr_ 



h.c. 



> (15) 



where Ck is a fc-dependent factor. Note that the term pro- 
portional to Ck only contains operators that couple states 
separated in energy by at least fko and result in correc- 
tions in eigenenergies oc fcp and can thus be dropped. 
This results in the effective Hamiltonian for the Rashba 



k 2 

H R ,cs = — 
1 



k R ka z 



41 



(2h+l)a z 



4 (1 + 2k R k) 



where the function 7 is defined as 
-ik R 



7 = 7(n) 



V2 



1 - 



4(1 + 2k R k) 



(16) 



(17) 



The order of the operators in the last term of Eq. 
important since [a, 7] 7^ 0. 
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B. Diagonalizing f?R,eff 

An interesting feature of Hr^s, Eq. (|16[) , is that it can 
be diagonalized via a generalized rotation matrix R(k) 
(see Appendix IA"| . which results in 



#R,diag — R H~R,eftR 

H+(n) 
H-{n) 



(18) 



where H± (h) are given in Appendix [X] The eigenstates 
of -f/R,diag(fc) are denoted by \k, n, s) and corresponding 
the eigenenergies £ns(fc) are given by 



4/2 



£ n,;(^) 



2 



1 



2A: R fc 



4/2 



1 + 2k K k 



A„(fc), (19) 
A n+1 (fc). (20) 



The generalized Jaynes-Cummings coupling mixes adja- 
cent transverse states with opposite spin, that cross at k* 
(indicated by the empty circles in Fig. [1} . This mixing 
is described by A„(fc) 



1 - 2k R k 



k^n 



+2k%n ( 1 



1 + 2k R k 

k R n 
4(1 + 2k R k) 



1/2 



Note that the (n, s) = (0, f) state reduces to 
k' 2 



— kvk 



2(1 + 2k R k) ' 



(21) 



(22) 



which reflects the fact that |0,f) is an eigenstate of the 
Hamiltonian in Eq. (|16l) . 

The above perturbative scheme does not depend on the 
fact that both W and V have equal coupling strengths. 
For completeness we give the eigenvalues for different 
coupling values in Appendix [Cl 



C. Numerical vs analytical results: Rashba case 

To better understand the range of validity of the 
above approximation scheme, it is instructive to look 
at the original Hamiltonian in Eq. When k = 

the Hamiltonian is easily diagonalized since it is just a 
shifted harmonic oscillator. The eigenstates of H R (0) are 
eT a y \n, s), here |0, n, s) — \n,s), with eigenvalues 



£n, s (0) = n - 



kx 



(23) 



Note that only the eigenvalues in Eq. (1221 reduce to the 
exact solution in Eq. those in Eqs. (HU) - (gUJ) do 

not. The magnitude of this deviation can be used to 
determine the accuracy of the analytical spectum. The 



origin of this deviation is that we have chosen to find an 
approximation scheme valid for k > 0, but not at k = 0. 
This allowed us to do the necessary approximations to 
obtain (I16[) , whose maximum deviation from the exact 
results occurs at k = 0. Comparing Eqs. (l20l) - ([22|) to Eq. 
()23[) gives the deviation of the approximate solution and 
the exact one at k = 0. Requiring that the deviation 
of the approximate solution and the exact solution be 
much smaller than the energy separation of the transverse 
states 



Ae = |e„.(0)-(n- *&/2)| < 1, 



(24) 



yields the inequality 1, where we have expanded 

the square root in Eq. (|21[) in powers of k R . Fixing the re- 
quired value of the absolute accuracy S, e.g., S = 5 x 1CT 2 
for a 5% accuracy, determines the number of transverse 
states that satisfy Eq. (1241) . The maximum value of n is 



< 



_ 2(46) 



1/3 



(25) 



which gives, e.g. n max = 18 (116) for fcrt = 0.25 (0.1) and 
5 = 0.05. In practice, this requirement is not very restric- 
tive since in most cases k R can be quite small and the 
approximate solutions are good for a very large number 
of transverse modes. More specifically, the values of a in 
usual III-V semiconductors range from 2 to 20 meV nm, 
which means that the above perturbation scheme applies 
to wires with £ < 1 /j,m and I < 100 nm, respectively. 

Since the spin-orbit coupling is usually a small effect 
the contribution of the kinetic energy fc 2 /2 is dominant 




FIG. 2. Eigenenergies of the first ten Rashba wire modes 
for kn = 0.25 (which yields k* =2). The solid lines are ob- 
tained within our analytical approximation scheme [Eqs. (|19p 
and (|20[) ]. while the empty circles correspond to the (exact) 
eigenenergies from a numerical diagonalization of Hn in Eq. 
0. The deviation between the numerical and the approxi- 
mate analytical spectra cannot be seen on scale of the plot. 
The dashed line denotes the Fermi energy Ef and the empty 
squares show the crossings of the states separated by 3&jq- 
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0.5 1 1.5 2 

k/k* 

FIG. 3. Similar to Fig. [2] but for the transverse states cor- 
responding to n = 50 — 59. Here the (Rashba) spin-orbit 
coupling strength &r = 0.10 =>• k* — 5. Note that even for 
such large n's the deviation of the analytical solution (lines) 
from the numerical (circles) can hardly be distinguished on 
the scale of the plot. 



turn heterostructures formed from semiconductors hav- 
ing bulk inversion asymmetry. The Dresselhaus spin- 
orbit coupling is present in III-V and II- VI material com- 
pounds, although its relative strength as compared to the 
Rashba coupling can vary. Using the same units of en- 
ergy and length as in Eq. Q the Rashba and Dresselhaus 
spin-orbit couplings (see Eq. ([!}) in our quantum wire 
can be written as 

#RD = ^ fc2 + atfl ~ (k-RVy ~ kv<r x ) k 

+-y=(a* - a)(kR(Tx - kucTy), (26) 

where ku is similar to k R with f3 instead of a. The pres- 
ence of the Dresselhaus term leads to extra terms which 
need to be taken care of before the proceedure of the 
previous section can be applied. In the next section the 
\a\ ^> |/3 1 case will be treated and in the subsequent sec- 
tion we will show how the \a\ -C |/3| and |a| rj |/3| cases 
can be written in the same form as used in the \a\ 3> \/3\ 
perturbation calculations. 

A. Couplings \a\ 3> 



in e ns (k) . To emphasize the details of the spin-orbit cou- 
pling we plot £ n s(k) — k 2 /2. The analytical solution of 
the spectrum (solid lines) and numerical solutions (cir- 
cles) for different values of the coupling strength are pre- 
sented in Figs. [2] and [3J In Fig. [2] the dispersion is plotted 
for k between and 3/c*. The dotted line shows the po- 
sition of the maximum value of the Fermi energy which 
we arbitrarily choose as Ep = 0.8(n max + 1/2), in order 
to insure that the error of approximation is less that 6. 
The crossing of opposite spin states separated by ~ 3Hujq 
starting at k « 3fc* are always far above Ep. This shows 
that the effect of discarding the third order coupling term 
does not affect the accuracy of the spectrum. The same 
is true for all other plots that follow and in those we 
restrict the value of k below 2k* . 

The transverse coupling mixes adjacent states with op- 
posite spin leading to a family of anticrossings, as seen 
in Fig. [5] This can be compared with Fig. [T] which shows 
the spectrum for k R = 0.25 and no transverse coupling. 
Note how the anticrossing at k* opens up for higher n 
leading to effectively flat bands for the low k part of the 
spectrum. This can be seen even more clearly in Fig. [3] 
for n = 50 — 59 and kn — 0.1, where the dispersion is 
very flat. This will result in the minima of e ns {k) moving 
from approximately k fcpj to k for higher values of 
n. 



III. RASHBA AND DRESSELHAUS 
COUPLINGS 

In addition to the Rashba spin-orbit coupling, there is 
also the Dresselhaus spin-orbit coupling present in quan- 



Let us introduce the angle cos 9 = k R / v^r + ^d an d 
the rotation operation 

kn<7 y -kncj x = yj fc| + fcp (<7 y cos 9 - a x sin 9) 

= k RB ei ea 'a y e-i 6cr \ (27) 

where &rd — v^R * s ^ ne effective spin-orbit cou- 

pling. Applying this spin rotation to the Hamiltonian in 
Eq. ([2"B)) , and then the rotation used in Eq. (jj) that takes 
{cr x ,cry) —> (cr z ,a y ), results in the Hamiltonian 

1 , , 

= —k + a' a — kjijyk a z 

+ !_^D( cos2 o . x _ S xn2ea z )(a f - a). (28) 
v2 

Note that the coupling becomes diagonal in <j z when 
cos 29 — 0. This occurs for values of the spin-orbit an- 
gle 9 = ±?, which corresponds to the case a — ±f3. In 
that case the exact spectrum reduces to shifted parabo- 
las. The term proportional to sin 29 can be removed via 
the unitary transformation 

U z = exp{-ik s xa z }, (29) 
where fc s = /crd sin 29 and x = ^=(a + a^). The result of 
applying the transformation to the Hamiltonian -Hrd is 
-Hrd = UIH^dUz 

= ifc 2 + a) a - k RD k<r z - 7j{k s ) 2 

(a+e-^V - a)e- tk °* - h.c. ) . 

(30) 
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Here k c — k^xi cos 29, along with k s , will play the role of 
perturbation parameters as k R did in Eq. ([5]) . The above 
Hamiltonian has a similar structure to that in Eq. ([5]), 
being identical to it when k s — 0, i.e., when 9 — 0,±5. 
The latter two values correspond to a pure Dresselhaus 
P = =p | ck | , where the sign is determined by the require- 
ment of the perturbation procedure: The velocity is pos- 
itive for s = — 1 states at k = 0, see Fig.[TJ The coupling 
in Eq. (|30p is written in terms of the exponential oper- 
ators in order to remove all o~ z contribution. Thus, the 
structure of the transverse coupling in the pure Rashba 
case is maintained, i.e. only adjacent states with opposite 
spins are coupled. By rewriting the coupling 



-ik s a 



l (a f ~a)e 



— ik s a 



a f e - lv / 2fc B a 
_ e -iv / 2fc s a t ae -iV2fc s O ) ( 31 ) 



and writing the exponential operators as a power series 
one can show that the diagonal matrix elements of the 
operator in the above equation is zero. The exponential 
operators add the complication that states with opposite 
spins separated by NHoj, N > 2, are coupled with cou- 
pling strength oc k^" 1 = fe^o , since we have assumed 
that sin 29 < £;rd- This assumption limits the possi- 
ble values of 9 but the range of allowed values of 9 can 
be extended by repeating above the arguments around 
\a\ « \P\ and \a\ -C \P\, which will be shown in the next 
section. 



1. Effective Rashba-Dresselhaus wire Hamiltonian 

Now we can proceed as was done with Eq. (jlip using 
the same form of S but using the substitution k R —> &rd 
in the o~ z term and k R — > k c in the transverse coupling, 



-s 



(32) 



Using the same arguments as following Eq. (|15|) we obtain 
the effective Rashba and Dresselhaus Hamiltonian for the 
wire by discarding all terms that lead to correction of 
order (fc RD ) 4 : 

„ r . (fc c ) 2 l + (2n+l)^ 

H RB , S = -+n- fc RD fc a z - 4 (1 - 2fcRpfc) 



+ ^(7RDa t c r + + acr_7^ D ) 



where 7jj/) is defined as 

-»(fcc) 2 

V2 



7rd 



1 



(k c ) 2 a*a 
4(1 + 2k RU k) 



(33) 



(34) 



Equations (|33|) and (|34| have the same form as the corre- 
sponding equations for the pure Rashba case, except for 
the (<7-fQ^ +(7_a 2 ) term. 



2. Diagonalizing HBx>,eB 
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0.5 1 1.5 2 

k/k* 

FIG. 4. Eigenenergies of the first ten transverse modes of 
the Rashba-Dresselhaus wire. Here the spin-orbit coupling 
strength fc RD = 0.25 and 9 = 0.126 (a > P) . The solid lines 
denote the approximate analytical solution in Eq. (|37[) . while 
the empty circles correspond to the numerical diagonalization 
of the corresponding Hamiltonian, Eq. (|26p . The inset zooms 
in on the n = 9 dispersion, showing more clearly the devi- 
ation between the analytical and numerical solutions. The 
maximum deviation is about 4% near k — 0. 



The goal now is to find a transformation that diago- 
nalizes Eq. (1331) . Instead of finding the exact transforma- 
tion we will take advantage of the fact that the rotation 
operator used in the pure Rashba case transforms the 




0.5 1 1.5 2 

k/k* 

FIG. 5. Similar to Fig. g] with fc RD = 0.10 and 9 = 0.126 
(a 3> P). The deviation of the analytical solution (lines) 
from the numerical (circles) is of the order of a few percents, 
and on the scale of the plot can hardly be distinguished. 
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n = 50 




k/k* 

FIG. 6. Similar to Fig.[5]but for the wire modes corresponding 
to n = 50 - 59. Here k RD = 0.10 and 6 = 0.126 (a > f3). 
The inset shows more clearly the deviation of the analytical 
solution (solid line) from the numerical one (circles) for the 
n — 59 mode. The maximum deviation is less than 4% near 
k = 0. 



squared ladder operator terms into an effective magnetic 
field term B(n)a x plus terms that can be treated per- 
turbatively. This proceedure is explained in Appendix 
iBl The diagonal form of the Rashba and Dresselhaus 
Hamiltonian is thus 

-ffRD;diag = 

_ t /ff+(n) B(n) . 
~ 1 B(h) H-(n) 1 



(35) 



where R is the same as in Eq. (fl8|) with renormalized 
parameters [Eq. (1551) ] and 1Z is the operator that diago- 
nalizes the 2x2 matrix (which commutes with h) in Eq. 
(|35| . The quantity B(h) is given by 



(36) 



B(n) = fc c fc s 2 ^ + 1) (l - cose(n)), 

see Appendices lAl and IB! and the eigenvalues of Eq. ([35 
are 



gryt + £rq ± . / (^n,t ~ £ n,\) 2 



B 2 (n). (37) 



Note the form is just what one would expect from adding 
a magnetic held along the x-axis to a system with a mag- 
netic field, and quantization axis along z. 



3. Numerical vs analytical results: Rashba- Dresselhaus case 

Figures [4] - [6] show the comparison of the analytical 
[Eq. (|3"T1) ] and numerical eigenenergies of the Rashba- 
Dresselhaus wire, for different values of /crd and the max- 
imum value of 9 = 0.126, reflecting that sin 29 < &rd- 



The numerical eigenvalues are obtained by diagonalizing 
the Rashba-Dresselhaus wire Hamiltonian in Eq. (|26p . 
The insets in Figs. 2] and |5] zoom in on transverse states 
n = 9 and n — 59, respectively. On this scale the dif- 
ference between the numerics (circles) and the analytical 
solution (solid line) is more evident. The maximum devi- 
ation, closest to k = 0, is of the order of 4% for the highest 
values of n but it should be noted that the relative error 
if further reduced by a factor w 59.5. Depending on the 
required accuracy, the analytical results may be extended 
beyond the strict requirement that sin 29 < fc^D ■ 



B. Couplings \a\ <C \(3\ and |q| ~ \/3\ 



The rotation introduced in Eq. (f28j) consists of two 
parts. The first part has the purpose of rotating the 
spin-orbit coupling term linear in k along a y , resulting 
in the same coupling as in the pure Rasbha case. The 



second rotation simply transforms a y 



<7, and docs 



not affect the perturbation scheme itself. The former 
rotation works for any value of 9, although the pertur- 
bation scheme for the transverse coupling only works for 
sin 29 < fcRD- With this in mind the cases |a| ~ \/3\ and 
| a | <SC |/3 1 can be treated in a similar way. 

When |a| <C |/3|, it is convenient to parametrize the 
spin-orbit angle as 9 — ^ + 5, where 6 is a small pa- 
rameter. The same steps are taken to arrive at the 
same equation as in (|2"51) apart from the new param- 
eters k s = fcRD sin(7r + 28) = — A;rd sin(2<5) and k c = 
/crd cos(7r + 25) = — /crd cos(2(5). Since the spectrum 
is even in the transverse coupling parameters fc CjS the 
eigenenergies for |a| <C |/3| are the same as those in Eqs. 

dSZD 

If the couplings are of similar strength, |a| ~ the 
angle is parametrized as 9 = ^ — 5, where 6 is again a 
small parameter. The Hamiltonian after the rotation is 



1 „ , 

-k + a 1 a - k RB k<7 z 

+ (sin 25a x — cos 28a z ) (a< — a) , (38) 
v2 



which is identical to Eq. (1281) apart from the switched 
roles of sine and cosine due to the phase difference of 
7r/2. As before, we can now proceed with proper unitary 
transformations in order to arrive at an effective Hamil- 
tonian, which can be similarly diagonalized to obtain an 
approximate analytical expression for the eigenvalues. 

Figure [7] shows numerical and analytical results for the 
eigenenergies of the Rashba-Dresselhaus wire with /crd = 
0.25 and different values of 9: 9 = ir/2—8 and 9 = ir/4—8, 
with 8 = 0.1. The solid line and the empty circles denote 
the 9 = tt/2 — 8 results of the analytical [Eq. (l37l) ] and 
numerical solutions, respectively. The difference between 
the analytical and the numerical results are similar to the 
case with 9 = 8, i.e. they are close to 5% corresponding 
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9.5 - 



K^sul^.. . o</? : (9,t) 




0.5 1 1.5 2 

k/k* 



FIG. 7. Here we plot the numerical and analytical solution 
for /crd = 0.25 with a <C ft (solid line and circles) and a ~ /3 
(dashed line and squares) for levels (5,t) and (9, f). The 
analytical solution is the same as that derived for the a <^ /3 
case. Note that the a « /3 solution is accurate near k = but 
deviates more for k ~ 2k* , see text for details. 



the estimate given at the end of Section|TTJ The 8 = ir/4— 
5 results look slightly different. The maximum deviation 
between the analytical and the numerical results occurs 
near k/k* — 2, where the contribution of the a 2 , a) 
terms is the greatest. 

This is best understood by looking at Eq. ([5]). The 
phase shift reverses definition of k s and k c so k c oc fcp D 
and fc s oc fcRj> Contributions from higher order terms in 
the exponential are thus larger and lead to the deviation 
around k/k* — 2. This demonstrates that the perturba- 
tion scheme for \a\ <C |/3| also works, within the accuracy 
determined by S, by adjusting the value of the spin-orbit 
angle according to |a| ~ |/3| and |a| ^> |/3|. 



IV. MAGNETOTRANSPORT IN RASHBA 
WIRES AND EFFECTIVE g FACTORS IN SO 
COUPLED 2DEGS 

In this section we treat two relevant physical problems 
in which the approximation scheme developed in the pre- 
ceding sections can be applied: (i) magnetotransport in a 
Rashba wire and (ii) the calculation of g-factors in spin- 
orbit coupled (Rashba-Dresselhaus) 2DEGs. In both (i) 
and (ii) we consider perpendicular magnetic fields. 

A. Shubnikov-de Haas oscillations 

In quantum wires the Shubnikov-de Haas oscillations 
in the magnetoresistance (along the wire) arise from the 
subsequent depopulation of the transverse wire modes as 
the bottom of their bands cross the Fermi energy Ep 




012345678 

k 



FIG. 8. B = eigenenergies of a Rashba wire [Eq. (|37[) ] as a 
function of k for a Rashba coupling strength fca = 0.135 and 
E F = 30&JO- 



for increasing magnetic fields 29 . This is similar to what 
happens in 2DEGs and also in bulk metal systems [see 
Chapter 11 in Ref. l30j . Similarly to 2DEGs, in quantum 
wires Shubnikov-de Haas oscillations allow the extraction 
of spin-orbit coupling strengths from the oscillation beat- 
ing patterns. Following the prescription in Refs. [3l| and 
[6|, below we outline in detail how magnetoresistance beat- 
ing patterns can emerge from the uneven distribution of 
the positions (in magnetic field) at which the subsequent 
bottoms of the transverse bands cross Ep. 

Figure [8] shows a plot of the dispersion relation [from 
Eqs. (fT9"|) and (|2TJ)) ] for a parabolic quantum wire with 
Rashba spin-orbit coupling /cr = 0.135 in the absence of 
a magnetic field. Note that here the plot includes the 
kinetic energy contribution. Applying a magnetic field 
B perpendicular to the 2DEG increases the effective con- 
finement of the wire, resulting in an effective confinement 
strength of 

n = ^l + ^ c , (39) 

where uj c = eB/m* is the cyclotron frequency. This 
raises the bottom of the bands, thus shifting the spectrum 
to higher energies. As the magnetic field increases the 
highest states below Ep are pushed above Ep. This gives 
rise to a peak in the magnetoresistance. As the magnetic 
field is further increased more states cross Ep thus lead- 
ing to oscillations in the Shubnikov-de Haas signal^ 2 -. As 
was discussed in Sec. [IT] the minima of s n . s (k) occur at 
k « 0, at least for high enough n. For lower values of n, 
where the minina is closer to k « k^, and the approxi- 
mation k leads to an error in the roots of order fc^. 
Hence, further below when we need to determine the val- 
ues of B for which the highest energy level crosses Ep we 
will use e ns (k = 0, B) = E F M± 

Using the ladder representation the Hamiltonian of a 
quantum wire in a perpendicular magnetic field at k = 



9 



(wave number along the wire) is given by 



H 



2V2 



r 2 a'a — 



ik 



4 mo 



R+ / + ^ 



(40) 



where r — uj c /ljq and the bare electron mass m appears 
through — = . The coupling coefficients are defined 

G mo 2e " 

as 



k R + = Wl + r 2 (l + r 2 /(l + r 2 )) (41) 

fc R - = Wl + r 2 (l - r 2 /(l + r 2 )). (42) 

In writing Eq. (|40l) we assumed a negative g-factor. By 
comparing Eqs. (f4"U| and (J5J) we see that g* plays the 
same role as k in the cr z term. The relation £ n ,s{k < 
0) = e n - s (k > 0) gives the spectrum for the opposite 
signs of the g-factor by simply reversing the spin label, 
as discussed following Eq. © on how to access the k < 
part of the spectrum. Note that the coupling strengths 
above are not identical, since the magnetic field breaks 
time reversal symmetry. This is in contrast to Eqs. ([7]) 
and ([5]). However, as we pointed out following Eq. © 
this does not affect the approximation scheme we have 
developed. The eigenvalues are obtained by substituting 
the parameters in Eq. (gOJ) into Eqs. (|C2|) and (|C3)) in 
Appendix ICl 

Figure[9]shows the bottom of the energy bands e ns (k = 
0, b) as a function of the magnetic field. As the magnetic 
field is increased the k = levels subsequently cross Ef, 
thus depopulating completely the corresponding trans- 
verse mode. These crossings occur at discrete values of r 
that we denote by fj, which correspond to magnetic field 
values Bi = ^f^n. The black dots in the inset of Fig. [3] 
shows the first 13 values of corresponding to solutions 
of £ ns (k — 0,ri) = Ef/luq for different values of (n, s). 
This results in a sequence of roots {fi}, each correspond- 
ing to the specific magnetic field value Bi — m e "° . 

Plotting the difference between subsequent roots bri = 
^"i-t- 1 — fi as a function of gives rise to the data (black 
dots) in Fig. [TO] The dashed envelope in Fig. \TU\ high- 
lights the beating pattern of the crossings. As dis- 
cussed in Ref. @ the data in Fig. [10] is directly related 
to the Shubnikov-de Haas oscillations of the wire: in- 
tuitively, the subsequent magnetic-field-induced depopu- 
lation of the transverse subbands affects the resistance 
of the wire because the number of conducting channels 
is being reduced. Hence the beating pattern describing 
the magnetic-field positions of the subsequent crossings 
of the band bottoms should manifest itself also in the 
magnetoresistance of the wire. 

We emphasize that the observed beating pattern arises 
from the superposition of the two spin-dependent contri- 
butions that have different frequencies: the Zeeman split- 
ting and the Rashba coupling. The parameters used in 
the figure are k R = 0.125, 0.135, and 0.145, g* = -4, and 
ro/mo = 0.037, corresponding to Gao.23Ino.77As. Assum- 
ing Hujq = lmeV, the value &r = 0.135 corresponds to 




« c /«o 



FIG. 9. Eigenenergies [Eqs. (fl9|) and pO)l ] as a function of the 
magnetic field with k = 0, ^rd = 0.135, g* = —4, m/mo = 
0.037, and Ef = SOHluq. The inset shows a zoom of the 
intersections of the dispersion curves with Ef- 




FIG. 10. Separation between subsequent roots Ar ; = r i+ i — 
Vi, plotted as a function of r^. The separation of the node 
points, e.g., r « 1.0 and r ~ 2.5 for the uppermost curve, is 
determined by the value a/(lohw). 



a Rashba coupling a = 6.38meVnm; both are realistic 
experimental values.—^ By examining the Shubnikov-de 
Haas oscillations obtained in experiments the difference 
in magnetic field between the node points can be ex- 
tracted. The value of the experimental node separation 
(in terms of magnetic field) can the used to extract the 
spin-orbit coupling from the numerics shown in Fig. 1101 
i.e., finding a that will give the node point separation 
equal to the experimental value. 

Since in our approach the determination of the relevant 
intersecting fields B t = m "" r-t is straightforward via the 
transcendental equation s ns (k = 0, r^) = Ep/^o), we 
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believe the procedure just outlined should make it sim- 
pler to extract of the Rashba spin-orbit coupling from 
Shubnikov-de Haas oscillations in quantum wires. 



B. Landau-level g- factor 

Our approximation scheme can also be used to obtain 
the effective g-factor of a 2DEG in a perpendicular mag- 
netic field and in the presence of both Rashba and Dres- 
selhaus interaction. Again assuming a negative g-factor, 
the Hamiltonian of the system is 

1 1 n 

H = ^(Ill + Il 2 v )--\g*\n B BG z 



2m 

4 

where II ■ 
operators 



h 



(U x a x — Tly&y 



< D 



r x e. 



(jly(T x — TL X <7y) 

(43) 

By defining the usual ladder 



V2H 



(u x + m y ), 



(44) 



where £ c = \J h/eB, the Hamiltonian can be written as 



H = huj c a} a 2~HbBct z 



20 



■(aV+ + aer_) 



i2a 
2V2£, 



(aV_ — aa + ). 



(45) 



Through the transformation a — > ae 



and <t+ 



CT+e '*'/ 4 and measuring the energy in terms of of the 
cyclotron energy huj c — heB/m* we recover the Hamil- 
tonian in Eq. ([5]) with renormalized parameters. Due 
to the magnetic field, the Hamiltonian in Eq. (|41d)) is no 
longer symmetric under time reversal and the Jaynes- 
Cummings term W and the perturbation term V are no 
longer transformed into one another under time reversal, 
but as mentioned before this does not affect the approx- 
imation scheme. 

The energy spectrum of the Hamiltonian in Eq. f|45[) is 
obtained by substituting the renormalized coupling pa- 
rameters into Eqs. (|C2l) - (|C4p . in AppendixICl This shows 
that the Hamiltonian for a 2DEG in a perpendicular mag- 
netic field in the presence of Rashba and Dresselhaus cou- 
plings is formally equivalent to a parabolically confined 
2DEG with Rashba interaction. Note that the spectrum 
for pure Dresselhaus coupling, a — 0, reduces to the 
known result for 2DEG in a perpendicular magnetic field 
and with Dresselhaus coupling 33 . 

From the energy spectrume ns the (/-factor can be de- 
fined for different Landau levels. We do not include the 
effects of electron-electron interaction, e.g. exchange en- 
hancement of the g-factor— , but rather focus on the effect 
of the spin-orbit coupling. The ^-factor of Landau level 
n is defined as the difference of Kramers doublets divided 
by hbB 



fjn 



fl B B 



(46) 



— I 1— 



- 1 1 — 



= 1, n = 3 




a [meV nm] 

FIG. 11. n t,l -Landau-level <? n -factors for two different values 
of \g*\ — 4 and \g*\ = 1, as a function of Rashba coupling. 
Other parameter values are B = 2 T, ^- = 0.037, and /? = 
6meV nm. 



In the absence of spin-orbit coupling the <?-factor is equal 
to g* for all Landau levels. 

In Fig. [TT] we plot the effective g-factor for the first 
four Landau levels for two values of \g*\. The parameter 
values are B = 2 T, ^ = 0.037, and (i = 6 meV nm and 
different g- factor values \g*\ — 4 and \g*\ — 1, while a 
ranges from zero to 8 meV nm. 

For the lower value \g*\ — 1 the effective <?-factor 
changes sign for n > 2 for low values of a. As a is 
raised all the g n 's cross at a single value of a. Similar 
crossing behavior is observed for \g*\ =4 but the sing 
of g n remains the same since for n = — 3. The results 
shown in Fig. [TT] are best understood by expanding Eq. 
(TTrJl) assuming that a, f3 <C huj c £ c /2. The ^-factor then 
becomes 



9n 

1.9* 



8(n + |) / (fc D 4) 5 



(fc R 4) 2 



m \ I — lull m I _i_ lull m 

2 mo V 2 mo 2 mo 



47) 



Note that <C 1 for typical parameters. The above 

equation shows two important things. First, the (7-factor 
can change sign when the spin-orbit part exceeds — \g*\ 
which can occur for high enough values of n or /3. Sec- 
ondly, the effect of the spin-orbit interaction is completely 
canceled when the Rashba and the Dresselhaus couplings 
satisfy the condition 



(fc D 4) 2 = (k R e c f 



1 



|g*l m, 

2 toq 



1 



|g*l to 

2 mo 



(48) 



Such crossings are a general feature of having both 
the Rashba and the Dresselhaus couplings but the sign 
change can only occur for small enough values of |g*|. 18 
Note that if g* = 0, this condition corresponds to a = 
±/3, where the energy spectrum is known and the spin 



11 



splitting of Kramers doublets is zero. The above results 
may prove useful for experimentalist to extract spin-orbit 
coupling strengths via g-factor measurements. 



V. CONCLUSION 

We have developed an analytical approximation 
scheme suitable for obtaining the energy spectrum of sys- 
tems with both the Rashba and Dresselhaus spin-orbit 
couplings and parabolic confinement, e.g., due to a mag- 
netic field or electrostatic gating. We have applied our 
approach to investigate (i) magnetotransport in quantum 
wires with the Rashba coupling and (ii) the effective g- 
factor of a 2DEG in a perpendicular magnetic field in 
the presence of both the Rashba and the Dresselhaus 
couplings. In (i) our approximate eigenvalues allow us to 
obtain beating patterns that relate to the Shubnikov-de 
Haas oscillations, by simply finding the roots of tran- 
scendental equations. For (ii) we derive a relation for the 
Landau level dependent effective g„-factor. We find that 
for small enough bulk g-factors the spin-orbit effects drop 
out of g n for certain values of the Rashba and Dresselhaus 
couplings. We believe our analytical approach should al- 
low for a more straightforward extraction of spin-orbit 
coupling strengths from Shubnikov-de Haas oscillations 
in wires and g-factor measurements in 2DEGs. 
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Appendix A: Rotation of generalized 
Jaynes-Cummings 

The Hamiltonians in Eqs. (JTCJ) and (|3"3"|) can be written 
in the form 



H, 



off 



__ / hi(n,k) j(h,k)a' 
aj*(n,k) h 2 (n,k) 



(Al) 



where h — a^a and the diagonal elements are defined as 
Til (0*0, k) = (+l|#eff| + 1) and h 2 {a)a, k) = (-l\H eS \ - 
1). All the matrix elements depend on k but since it is 
only a parameter it is dropped for convenience. The goal 
is to find a unitary transformation, or rotation, R that 
diagonlizes the Hamiltonian in Eq. (|A1[) . If the ladder 
operator a((r) is replaced by a complex number a (a*) 
the rotation operator reduces to 



Ro 



hi 



hi 



7=== sin 

/a* a ^ 



7=5= sin S 

/a a £ 
COS % 



(A2) 



Note that the factor . . 7 , a can equivalently be written 

\~f\\/a*a ^ J 

in the ususal phase factor form e 1 ^. The choice of our 
notation will help to understand the form of the general 
rotation matrix discussed below, see Eq. (|A5|) . To take 
into account the commutation properties of the ladder 
operators in the diagonalization of Eq. (|A1[) . we start 
with the form 



R = 



-X 



y/l+XXf 



(A3) 



y/l+XXi y/l+XX* 



where X is an operator that depends only on the ladder 
operator, to be determined. The diagonalization requires 
that the 2x2 spin matrix 



Hd — R HefiR 



(A4) 



has zero off-diagonal elements. Using the ansatz X^ = 
— a^x(^) an d the relation f(n)a — af(h — 1) the diago- 



nalization is achieved when x(n) = —, 
the rotation operators takes the form 



: tan 



Hi,. 



and 



R 



(1 - Qo) + cos 



e(w-i) 



, 7*(") 

'h(«)l 



0(r, 



(A5) 



The diagonalization is achieved when 

e(A) 



ator Qq 



tan \ ' . Here we have introduced the oper- 



as —r 



a+a+1 



a. The operator Qq acts as the 



complement of the projector into the oscillator ground- 
state V = |0)(0|, i.e. Q \n) = |n) if n > and Q |0) = 0. 
It does not appear in any of the other matrix elements 
since only the operator 1/V1 + X^X gives rise to Qo, 
not the other form 1 / vT-j- XX*. This asymmetry comes 
from the fact that |0, +1) is an eigenstate of Eq. (|A1[) 
but not of 1 0,-1) or of any other state |n > 0, s). The 
angle operator 9(n) is defined via the ususal relation 



cos 8(n) 



h 1 (h + 1) - h,2(n) 



V^M" + 1) - h 2 (n)) 2 + 4(n + l)\-f(h + 1)| 2 

(A6) 



The form of the matrix in Eq. (|A5j) is reminiscent of the 
form of Eq. (|A2|) . the difference coming from the non- 
commutivity of a and a' and the operator Qq. 

The eigenenergies are determined by the diagonal ele- 
ments of iJd, which are given by 

[Hd] 2 , 2 = l((hi(h + l)-h 2 (h)) 

-y/faih + 1) - h 2 (h)Y + 4(n + l)| 7 (n + 1)| 2 ) 

(A7) 

N M = &i(ft)(l - Go) + - h 2 (n - 1)) 

+ ^(h 1 (h) - h 2 (h - l)) 2 + 4fi| 7 (fi)| 2 ) Qo- 

(A8) 
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The two equations above are obtained by algebraic ma- 
nipulations of Eq. (IA4|) and collecting all powers of a 
from l/Vl + X<X and writing them in terms of Qo = 



a* i 1 1 -i a and 1 



Appendix C: The eigenenergies for general coupling 
values 



For completeness here we provide the eigenvalues of a 
generic Hamiltonian of the form 



Appendix B: Rotation of o+a) + h.c. 



The effective Hamiltonian in Eq. (1551) cannot be ex- 
actly diagonalized using the rotation in App. [S] This is 



due to the term i (cr_| 



h.c. ). From the Hamiltonian 



that results from the rotation it is easy to discard terms 
which give corrections ©((/crd) 4 )- Using the property 
7 rSi = i, the transformed second order coupling is 

R^{a+a^ 2 + a^a 2 )R 

r— . 6(n-l) . 0(n) 

= \/n{n + 1) sin sin — — o x 

+ . 0(n- 1) 6(n) , 

+icr z I a 1 sin cos h.c. 

2 2 



1 \ ( t e (")\ , 
- I cr+ la 1 cos— ^— ) + h.c 



(Bl) 



As can be seen in Figs. [2] and |3j the only states that 
cross are the same transverse states with opposite spin, 
e.g. when n = 1 the f and 4- states cross at k ps 3.75 in 
Fig. [2] Taking advantage of this behavior we can deter- 
mine which terms in Eq. (|B1|) need to be kept and which 
can be dropped. The first term in Eq. (|B1[) couples states 
with with opposite spin but same n and the latter two 
terms couple adjacent transverse states that do not cross 
and can be treated perturbatively. They give rise to cor- 
rections in energy ©((/crd) 4 ) and can be discarded. From 
these considerations the form of the rotated second-order 
coupling is 



(cr+a 1 ' + er_a 2 )i? 



(B2) 



fccfc , Vn(n+l) (1 _ COS0( - )K + 



H 



+ a) a — Xa z + i 



Y 

2V2 



(aV 4 



(CI) 



+i — =(a^cr_ — &er+), 



where X, Y and Z are real numbers and the last term, 
proportional to Z, is the perturbative coupling. The only 
restriction on the parameter values are that X > and 
Y and Z should be of the same same order of magnitude 
although their numerical value can be different. This 
type of Hamiltonian covers all models discussed in this 
paper. Up to corrections of the order Z 4 , the resulting 
eigenvalues are 



, . v k 2 Z 2 /2 
£ n,tW = y + n- IW + A 7 



(C2) 



k 2 Z 2 12 

^n,x{k) = — +n+l-—^-^ n+1 , (C3) 



where 




+2Y 2 n 1 - 



, x 2\ 1/2 



Z l n 



4(1 + 2X), 



(C4) 



When Z = the spectrum reduces to the known Jaynes- 
Cummings solution. 
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